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1. Introduction 

The study of pseudospectra of operators is a very active area of research, of 
importance in applied mathematics and numerical analysis. We recall that if 
e > 0, the e-pseudospectrum of an operator Q acting on a Hilbert space is the 
set of complex numbers, A 6 C, such that 



|| (Q - XI)' 



> 
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or, equivalently, the set of A such that 



The importance of this set in applied settings stems from the following two facts: 

(1) If e is very small it may be difficult to distinguish the e-pseudospectrum 
from the spectrum of Q, 

(2) If Q is strongly non normal the e-pseudospectrum generally is much bigger 
than the spectrum, even if e is very small. 

In case Q is a differential or, more generally, a pseudo differential operator with 
a small parameter, h (such as a Schrodinger operator), it is natural to consider 
the asymptotic behavior of the e-pseudospectrum of Q where e is related to h, e.g. 
e = 0(h°°). Dencker, Sjostrand and Zworski have recently studied this problem by 
microlocal techniques, |}JJ- Earlier results for Schrodinger operators were obtained 
by E. B. Davies, jg|, and P. Redparth, QS]. 

On the other hand, in applications one often deals with large matrices, and then 

it is natural to estimate the pseudospectrum in terms of the size of the matrix. In 

a recent paper, ^7j; Trefethen and Chapman considered this problem for matrices 
T (N) = ( T Af) where 

(1-2) TjP = f { i- j)m odN(J/N), 1<3,1<N. 

Here the fj are 1-periodic coefficient functions . The main result of Trefethen 
and Chapman is that, under certain assumptions including the following "twist" 
condition, 

(1.3) <° 

where f(x,p) = ■ fj{p)e^ x , then A = f(xo,po) is in the e-pseudospectrum of Tjv 
where e = 0(e~ cN ) for some c > 0. 

The purpose of this note is to show that microlocal techniques can also be applied 
to the study of the pseudospectra of matrices such as (|1.2|) (and generalizations). In 
this light we interpret the twist condition as Hormander's solvability condition 



(1.4) {3?/, Qf}(x ,p ) <0 

on the Poisson bracket of the real and imaginary parts of the symbol of a pseudo- 
differential operator. Indeed if we define the Poisson bracket of the variables x and 
p above to be one, it is easy to check that is precisely Hormander's condition. 
The connection between Hormander's condition and pseudospectra was first made 
by M. Zworski in [2D] . 

In this paper we construct pseudomodes for Berezin-Toeplitz operators under 
condition (|l,4j) on the (smooth) symbol. Our construction is symbolic: in §2 
we introduce spaces of Hermite distributions containing the pseudomodes. From 
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the point of view of the symbolic calculus of these distributions, condition Q1.4[) 
is exactly the condition on e for the operator: 4- + ex to have a kernel in the 
Schwartz space of K., namely e > 0. 

Although we will discuss our results in detail in the next section, we should men- 
tion some limitations of our work. The methods of Trefethen and Chapman apply 
to rough symbols, /, and they obtain exponentially small error terms. For ana- 
lytic symbols, it is very likely that exponentially small estimates (in the Toeplitz 
setting) can be achieved by microlocal methods, as has been done in |Hj for pseudo- 
differential operators. The problem of dealing with general non-smooth symbols, 
however, is much more challenging. Trefethen and Chapman's main theorem in- 
cludes a global condition on the symbol (in addition to I1.2JI . and they present 
compelling numerical evidence that global conditions on non-smooth symbols are 
necessary for the existence of "good" pseudomodes (see §8 of |IZj). This is a very 
interesting issue that we do not address here. On the other hand, our results for 
smooth symbols are fairly general and include a number of cases not covered by the 
results in (e. g. the "Scottish flag" matrix). Furthermore, the pseudomodes we 
construct are localized in phase space, sharpening the localization results of |17j . 

We also mention that more straightforward microlocal methods can be applied 
to the study of non-periodic versions of (|1.2ft . along the lines of the example in 
§4.1. More generally, the Berezin- Toeplitz operator calculus opens up the entire 
spectrum of phase-space methods to study other problems associated with certain 
sequences of large matrices, and we hope to return to some of these problems in 
the future. 

1.1. The main results. The general setting for B-T operators is a Kahler man- 
ifold, X, together with a holomorphic hermitian line bundle, L — > X whose cur- 
vature is the symplectic form onl. If / : X — > C is a "classical Hamiltonian" , (a 
smooth function) consider the sequence of operators T f = {T ( f N) , N = 1, 2, . . .}, 

acting on the space Hn of holomorphic sections of the tensor power L® N , defined 
by: 

H N 3^^U N (f^), 

where ITv : L 2 (X, L® N ) —>TLn is orthogonal projection. The sequence, Tf, is the 
primary example of a Berezin- Toeplitz operator. More generally, one can allow / to 
depend on ./V as well, provided the TV-dependence admits an asymptotic expansion 
as N — > oo 

oo 

f(x,N)~J2 N ~ j M x ) 
j=o 

in the C°° topology. The function /o is then called the principal symbol of the 
operator. We include explicit examples of all this in §4. General recent references 
for the theory of Kahler quantization and Berezin- Toeplitz operators are [3], [7], 
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For each N Tljsr is finite-dimensional, and for N large, by the Riemann-Roch 
theorem, dimTij^ is a polynomial in N of degree one-half the dimension of X, 
n := ^dim X. 

The simplest cases are when X is a either the torus or the complex projective 
line, for which n = 1. Thus the parameter, N, is essentially the dimension oi TCn. 
Moreover, for such X the spaces TLn have a natural multiplicity-free representation 
of the circle group, whose eigenvectors form a canonical basis of TCn- We will write 
down explicitly the matrices corresponding to a B-T operator on these spaces in 
the next section. It turns out, for example, that the matrices (jl.2[) are the matrices 
of B-T operators on X equal to the two-torus. 

Our main result is: 

Theorem 1.1. Let Tt = {T^ N ' , N = 1, 2, . . .} be a Berezin-Toeplitz operator with 
smooth principal symbol f : X — ► C. 
0. For all AsC, 

\\(t\ n) -a/)(V0II r- 

inf f = inf |/(x) - A + 0{1/VN). 

■4>&n N \\ib\\ xgx 



1. Assume that A = / (xo) where xq £ X is such that 
(1-5) {K/,9/}(x )<0. 

Then there exists a sequence of vectors {i^n £ T~(-n} with microsupport precisely 
{xq} and such that 

\\(T ( f N) -AI)(^v)|| 
(1.6) f „ , ' = 0{N-°°). 

\\Vn\\ 

2. On the other hand, if A = /(xo) and 

(1-7) {K/,9/}(x )>0, 

then any sequence {V'Ar £ Wat} smc/i that M.b]) holds has microsupport away from 

We say a few words about the definition of microsupport in H2.ll 
In the case of more than one degree of freedom (i. e. if the dimension of X is 
greater than two), there are multiple pseudomodes under condition Q1.5JI . The 
level set / _1 (A) is a symplectic manifold (at least near xo) and we construct pseu- 
domodes associated with any germ of isotropic submanifold of it containing xo- 

We will also prove an additional result, analogous to Theorem 4 in whose 
hypotheses hold typically in case A is on the boundary of the image of the principal 
symbol. Let 5R(/) = /i, 9 ; (/) = f2, and for / C {1, 2} m denote by fi the repeated 
Poisson bracket: 
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where H g denotes the Hamilton vector field of g. Denoting the order of the Poisson 
bracket by |/| = m, we define the order of a point x G X as 

(1.8) k(x) := max{j G Z ; / 7 (sc) = for all |/| < j }. 

Theorem 1.2. Let A G dlmage(f) be such that: 

(1) dfx / for every x G 

(2) The maximum order k := maxj-gj-i^ is finite. 
Then there exist C, C\ > suc/i i/tai 

(1.9) dN- 1 ' 2 ^ inf J > CN'm. 



N 



The first inequality in (J1.9|) follows from Part of Theorem 11.11 the second one 
follows from subelliptic estimates. In the course of the proof of this Theorem we 
will also show that, in general, the microsupport of a sequence of vectors ipN £ 
TLn minimizing the Rayleigh quotient in ()1.9|) (i.e. the optimal pseudomodes) is 
contained in /~ 1 (A), see Proposition 13.31 

The existence part of Theorem 11.11 is proved by constructing pseudomodes out 
of a class of distributions that will be defined in the next section. The proof of 
Theorems 11.11 and 11.21 appear in §3, and §4 is devoted to examples. We present 
some additional results, including a description of the limit of the numerical range, 
in §5. 

Acknowledgments: We wish to thank Nick Trefethen for sharing with us an 
early version of and for encouraging remarks. Many thanks also to Maciej 
Zworski and to Thierry Paul for helpful conversations. 



2. Preliminaries 

2.1. Setup and strategy. Let L — > X be as in the previous section, and let 
P C L* denote the unit circle bundle in the dual of the line bundle L. We denote 
by 

(2.1) L\P) = ^Ll(P) 

the Fourier decomposition of functions on P under the action of the circle; explicitly 
/ G L 2 (P) is in L 2 k {P) iff f(e ls ■ p) = e tks f(p). We will also need the spaces 

(2.2) C?(P)=L 2 (P)nC»°(P). 

It is a tautology that C^°(P) (resp. L 2 (P)) can be naturally identified with the 
space of sections C°°(X, L® k ) (resp. L 2 (X,L® k ). We will henceforth identify these 
spaces without further comment. 

P is a strictly pseudoconvex domain, and under the natural action of the cir- 
cle group its Hardy space, TC, splits into Fourier components that are naturally 
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isomorphic to the spaces of holomorphic sections Ti n '■ 

oo 

(2.3) n =@H N , H N = H°(X,L® N ). 

N=0 

We denote by n : L 2 (P) -> H the Szego projector of P, and by ILv : L 2 (P) H N 
the orthogonal projection onto the summand TIn- 

The precise structure of the singularities of II has been known for some time, 
thanks to work of Boutet de Monvel and Sjostrand. We now recall the microlocal 
structure of II, as described in [I]. Let a denote the connection form on P, and 
let Z cT*Pbe the manifold 

Z = { (p, ra p ) ; p£ P , r > }. 

This is a conic symplectic submanifold of T*P, and II is a Fourier integral operator 
of Hermite type associated with the canonical relation 

Z A : |iU!:0 Z}. 

We will say a few words below about the symbol of II, referring to [I] for the general 
theory of Fourier integral operators of Hermite type. (See jTH] for a description of 
II as a Fourier integral operator with complex phase.) 

The overall strategy of our proofs is the observation that much of the asymptotic 
behavior of a sequence {ipx G C°°(X, L® N )} is encoded by the singularities of the 
distribution on P, tp = Y1n=i ^ ^ C~°°(P). For example, we have the following 
elementary result: 

Lemma 2.1. Given a sequence of vectors ipx G C°°(X, L® ), letifi = ^a?=i ^JV G 
C~°°(P). For each set, let H^ s -j(P) denote the Sobolev space on P consisting 
of distributions, u, such that (Ap + I) s / 2 {u) G L 2 {P), with the norm ||lt||( a ) = 
||(Ap + I) s I 2 (u)\\l2 . Then the following are equivalent: 

(a) i> G C°°(P) 

(b) For allseR Un\\( s ) = 0(i\r°°). 

If, in addition, ipw G Hn for all N, then the above are equivalent to: 

(c) H^vll^ = 0(N~°°). 

Proof. The spaces C£°(P) are invariant under Ap and orthogonal in Hr s \, so for 
each s = SatH^nH^)- It follows that (b) implies that ip G H^(P) for all 

s, and therefore it implies (a) (and (c), of course). Assuming (a) now, consider 
Dgip where k is a positive integer. Since this function is smooth it is in H^(P) 
for each s, and therefore = J2 N N 2k HV'ivll?^ < °o> which implies that 

||V?jv||( a ) = 0(N- k ). Therefore (a) implies (b). 
Let us define an operator A^ by the identity 

(2.4) Ap = A h + Dq. 
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Then [A^, Dg] = and the restriction of A^ to C^(P) agrees exactly with the 
Laplacian on C°°(X, L® N ) associated with the connection and the Hermitian struc- 
ture on L® N . This has the following consequence: If ipN S Hjv C L 2 (P), 

(2.5) Un\U = (iV 2 + iV + l) s/2 ||V7v||o. 

Indeed elements in TCn are eigenfunctions of Ap: We claim that Ap(ip]y) = N(N+ 
l)tpN- By virtue of ()2.4|) . this statement is equivalent to: A/^tv = Nipjy. This 
follows from the well-known Bochner-Kodaira relationship the metric and the d 
Laplacian on sections of L N , see for instance |1U] Proposition 6.1. Clearly (|2.5j) 
has the consequence that (c) implies (b) for sequences of vectors in TCn- D 

We end this subsection with a reminder of the notion of microsupport in Kahler 
quantization. Let {ipN G TCn} be a sequence of holomorphic sections of the tensor 
powers of L, and let tp = J^n^n- Since U(ip) = tp, the wave-front set of ip is 
included in Z. The microsupport of the sequence is defined as the subset of X 
which is the projection of WF(^): we say that x G X is in the microsupport of 
{ipN } if and only if 

3p x G P such that ir(p x ) = x and (p, a p ) G WF(-0). 

It follows from the above that microsupport of the sequence is the empty set 
iff H^aHIl 2 = 0(N~°°). In addition, one can show that xo G X is not in the 
microsupport of {i/jn G TCn} iff there exists a neighborhood, V, of xq such that 
sup xg y \iI>n(Px)\ = 0(N~°°), where for all x p x G P denotes any point projecting 
to x. 

The microsupport has a characterization in terms of the action of Toeplitz oper- 
ators analogous to the characterization of the ordinary wave- front set by the action 
of pseudodifferential operators. We refer to jjj, §5, for alternative descriptions of 
the microsupport. 

2.2. Polarized Hermite distributions. In this section we define and analyze 
the concept of generalized wave packets and their symbols in the context of Kahler 
quantization. In fact we'll define more general states, associated to isotropic sub- 
manifolds of a quantized Kahler manifold (although everything we do generalizes to 
any almost Kahler manifold quantized by a projector, II, with the same microlocal 
structure as the Szego projector.) 

Definition of polarized Hermite distributions. 

Let us begin by considering a closed conic isotropic submanifold 

KCZ. 

Obviously TZ is isotropic in T*P \ {0}, and therefore associated with it are spaces 
I l (P,1Z) of Hermite distributions on P. The general theory of such distribu- 
tions (together with many applications) was developed by Boutet de Monvel and 
Guillemin, see (see jl], or the Appendix for additional remarks). The polarized 
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Hermite distributions associated with 1Z are simply the projections of elements of 
I l (P,7Z) by the Szego projector: 

Definition 2.2. The space of polarized Hermite distributions of order I associated 
with 1Z is 

(2.6) I l n (P,K):=Il(I l (P,K)). 

We should point out that the composition Theorem 9.4 of [1] one has the inclu- 
sion: i l n {P,n) c i l (P,n). 

Notice that one has a natural isomorphism 

Pxl+ ^ Z 
(p, r) i ► (p, ra p ) 

which becomes a symplectomorphism if we put on P x IR + the symplectic structure 

(2.7) —d(ra) = —rda — a A dr 

(here r is the coordinate on the M + factor). Using this description of Z it is easy 
to show that the base of the cone, 1Z, is a submanifold Y C P such that the 
infinitesimal generator of the circle action on P is never tangent to Y . Moreover, 
the projection, it : P —* X, restricts to an isotropic immersion of F in X. Let us 
denote by Y this immersed isotropic submanifold of X. 

Definition 2.3. A sequence of holomorphic sections, {ipN G 7~Cn}, of the tensor 
powers of L will be called an Hermite state associated with Y iff there exists 
ijj £ ^n(-P) Tty such that the Fourier components of tp according to ()2.3|) are precisely 
the tpN. 

For example, if Y is a single point then the coherent states at that point are 
an Hermite state. The case when Y is Lagrangian (and hence Y Legendrian) was 
considered in 0. 

Symbolic matters. 

Our next step is to define the symbol of a polarized Hermite distribution. We 
begin by recalling the nature of the symbol of general elements in r{P,TZ). Such 
distributions have symbols which are symplectic spinors associated with 1Z. The 
definitions are made in the tangent space to T*P, so for notational convenience 
for each p G 1Z we'll set 

R p = TplZ, Z p = TpZ. 

Let 

N p := Rp/Rp 

where R° denotes the symplectic orthogonal of R p inside T p (T*P). N p is a sym- 
plectic vector space, called the symplectic normal space to 1Z at p. Abstractly, the 
symbol of u G (P, 1Z) at p is a smooth vector in the metaplectic representation of 
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the metaplectic group of N p tensored with a half-density along R p , i.e. an element 
of: 

Spm(R p ) := A l/2 (R P ) ® H^N,). 

Lemma 2.4. Let E p = {v £ Z p ; \/u £ R p uj(u,v) = 0}/R p be the symplectic 
normal of 1Z in Z at p. Then, the symplectic normal, N p , is naturally isomorphic 
to the direct sum 

(2.8) N p = E p © Z° p 

where Z° is the symplectic orthogonal of Z p inT p {T*P). 
Proof. This follows from the fact that 

r° p = w*®z; 

where R° z is the symplectic orthogonal of R p inside Z p , which itself follows from 
the fact that Z is a symplectic submanifold of T*P. □ 

It is important to note that the projection T p (T*P) — > T W ^X induces a sym- 
plectic isomorphism Z° ~ T n ^X. On the other hand, because of the negative 
sign in (|2.7|) . the same projection takes E p to (T^^Y /T^^Y)" , where the minus 
indicates a reversal of the symplectic structure. 

It follows from Lemma f2 . 41 that the metaplectic representation of the metaplectic 
group of N p is a tensor product: 

H(N p ) = H(E p )®H(Z° p ) 

(Hilbert space tensor product). It is this decomposition that reveals the structure 
of the symbol of a polarized Hermite distribution. In order to discuss this structure, 
we recall that the smooth-vector factor of the symbol of II is of the form e <g> e, 
where e G H^Z") is a normalized "ground state" (which can be identified with 
the ground state of the harmonic oscillator on T^i p \X defined by the metric). 

Proposition 2.5. The symbol of a polarized Hermite distribution, u £ Iji{P,1Z), 
is of the form: 

o- u = v u <8> k u <8> e, 
where e E H^Z®) is the symbol of the polarization and 

v u e f\ 1/2 (R p ) k u e H^Ep). 
By dividing a u by e one obtains the non-trivial map in the following exact sequence: 
_> f-Wftn) I l n (P,K) -» A 1/2 (R P ) HooiEp) 0. 
We relegate the proof of this technical proposition to an appendix. 
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3. Proofs 

3.1. Proof of Theorem EH Let T f = {T^, N = 1,2,...} be a Berezin- 
Toeplitz operator with smooth principal symbol / : X — > C. 

Part (0) of Theorem II. H is not difficult. Without loss of generality we can assume 
that A = 0. Since {(T N )*T( N ^} is a Toeplitz operator with symbol \f\ 2 , one has 
that for all tpN £ TLn, 

\\t^ n ^n\\ 2 = ((T N yT( N ^ N ,^ N ) ~ (u N (\f\^ N ),^ N ) = (\f\ 2 i> N ,M- 

More precisely, from the definition of B-T operators we have that for any sequence 
of ip N £ Hn, 

I|T (7V VaH| 2 = / \f(x)\ 2 \^ N (x)\ 2 d^ x + U N \\ 2 -0(l/N), 
Jx 

where d\i is the measure on X and |V>A"(aO| 2 is the square of the length of i^n{x) 
in the Hermitian norm of L® N — ► X. It follows that 

\\M N) - xi)(ij N )\\ r - 

inf f - ■ ,, > inf \f(x)\ +0(1/VN). 

*p N en N \{o} \\ip N \\ xgx 

In the other direction, let xq be the point of X where inf is attained, and 

let ifp^ = H(-,po) be a coherent state at a point po € P that projects down to xq. 
Then 

i> N en N \{0} \\WN\r WpoW 

We now prove part (1) of the Theorem. We begin with a few preliminary 
considerations. By Proposition 2.13 of [Tj. there exists a classical pseudodifferential 
operator of order zero on the circle bundle P, Q, such that: 

(1) Q commutes with the Szego projector ([II , Q] = 0) and with the S 1 action. 

(2) For each TV" the restriction of Q to Hn , Q '■ T~(-N —> Hjv is equal to T^ N \ 

(3) The principal symbol of Q satisfies: 

V(p, ra p ) e S <tq(p, ra p ) = f(ir(p)) 
where ir : P — ► X is the projection. 
Assume now that A = /(xo)> where xo G X is such that 
(3.1) {K/,3f}(x )<0. 

The inverse image / _1 (A) is, in a neighborhood of xq, a codimension-two symplectic 
submanifold of X. Let us pick an isotropic submanifold of / _1 (A) (not necessarily 
closed), Y, containing xq. Then Y is an isotropic submanifold of X. Our consider- 
ations are local: we restrict our attention to a neighborhood of xo where {3i/ , 3°/} 
is negative and such that there exists a lift of Y to a conic isotropic submanifold, 
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1Z C Z, in the sense of the previous section. We will construct a pseudomode with 
microsupport equal to Y . Notice that we may take Y = {xo} if we wish. 

Since [II, Q] = and II is self-adjoint, we have [II, QT] = o. Therefore, by 
Proposition 11.4 of 0], 

(3.2) V(p,ra p )GS {Ka Q , %a Q }(p, ra p ) = {&/ , 3?/}(tt( P )), 

where oq is the principal symbol of Q. Notice that the Poisson bracket on the left 
is on T*P (with respect to the cotangent bundle structure), while the one on the 
right is the Poisson bracket on X. Therefore, the Poisson bracket conditions (|3.H 
13 .5|) on / are inherited by <tq. 

To see how the Poisson bracket condition becomes relevant in our considerations, 
we first consider a calculation in the Heisenberg representation on L 2 (M. k ). 

Lemma 3.1. Let C be the operator on L 2 (M. k ) corresponding to the action of 
£ £ (M. 2k ,u>) (g> C under the Heisenberg representation. Then if 

>o, 

then the restriction of C to the smooth vectors maps S(JH k ) onto itself, with a 
non-zero kernel. 

Proof. The metaplectic representation describes how C transforms under the action 
of the symplectic group on £. That is, C g ^ = U (g)C^U (g -1 ), where g i— > U(g) is the 
projective unitary representation that gives rise to the metaplectic representation 
when we take the double cover. The action of the metaplectic group preserves the 
smooth vectors 5(M fc ), so in our argument we can replace £ by g.£, for g symplectic. 

Under the assumption that e := u>($i!;,Qt;) > 0, it is a straightforward exercise 
to see that g can be chosen so that = ee\ + if%, where ei, . . . , e&, /i, . . . , is 
the standard symplectic basis for M. 2k . Thus it suffices to prove the result for 

L = h exi. 

ox i 

We see then that ker £ n 5(M fc ) contains functions of the form ip(xi, . . . , xp.) = 
e - ex l/ 2 a{x 2 , ■ ■ -,x k ). And to show that C maps S(R k ) onto itself, let / G S(R k ). 
The ODE Cu = f can be solved by variation of parameters: 

u( Xl ,...,x k ) = / f(t,x 2 ,...,x k )e e{t2 - x " )/2 dt. 
Jo 

For e > 0, the estimates showing that u E 5(M fc ) follow easily. □ 

Applied at the symbol level, Lemma 13 . 1 1 leads directly to the following construc- 
tion: 

Proposition 3.2. There exists a distribution u in the class In(P, Tl) of polarized 
Hermite distributions of order zero associated with 1Z such that: 



(3.3) 



(Q-xi)u e c°°(P). 
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Proof. Suppose u G I^(P,TZ) with a u = v u ® k u ® e as in Proposition ^. 51 It's clear 
that {Q — XI)u = n(<5 — XI)u is also a polarized Hermite distribution. What is its 
symbol? Since the principal symbol of Q — XI vanishes on 1Z, we are led to use the 
first transport equation for the Hermite calculus. The Hamilton vector field, £, of 
the symbol of Q at a point p G 1Z is in the symplectic normal space N p . Therefore, 
the Heisenberg representation of that space associates to £ p an operator, £, on 
the space H^Np). According to Theorem 10.2 in @], (Q - XI) (u) G ^^(P^U) 
and its symbol is v u <g> C(k u (g> e). In fact, since (Q — XI)u is still polarized, under 
the decomposition H OQ (N p ) ~ H OQ (E p ) H^Z ), L acts only on H^iEp). The 
symbol is really v u ® C(k u ) (g) e. 

We noted in the previous section that the pull-back of the symplectic form from 
X under the natural projection E p — > T^f p \X is the opposite of the symplectic 
form on E. So the Poisson bracket condition (j3.2j) along with Lemma 111 II implies 
that, as an operator on H OQ (E p ), C is onto and has a non-trivial kernel. 

Because C has a kernel, we can choose uq G I°(P,1Z) with symbol C(a(uo)) = 0. 
Therefore v\ := (Q — XI)(u\) G I~ 1 {P,1Z). Because C maps onto H^Ep) (by 
Lemma Hill again), we can then find U\ G I^ 1 ^ 2 {P,TZ) such that C(a Ul ) = —o~ Vl , 
Thus V2 := (Q — XI) (uo + u\) G I~ 3 / 2 (P,1Z). Continuing in this fashion and 
finishing with a Borel summation of the Uj's produces u G I°(P,H) such that 
(Q — XI) (u) is of order (— oo), and therefore smooth. □ 



We can now finish the proof of the part (1) of Theorem 11.11 Let u be as in 
the previous Proposition, which we choose to have a Gaussian principal symbol. 
Define um '■= Hjsr(u). By Lemma 12. 11 (j3,3|) implies the norm estimates 

(3.4) ||(TW - XI)u N \\ = 0(N-°°). 

On the other hand, letting tpjfp be the coherent state at po G P, from the repro- 
ducing property un(po) = (tPpP,UN) we obtain the estimate 

ii I, ^ Wn(po)\ 

W U n\\ > (AT),, • 

Wpo II 

Combining 1)A.1|) with the asymptotics 

\\<Ppo II =njv(po,po) ~ (^J 

(see e. g. equation (28) in &\), we see that 

IMI > cn-^'\ 

where I = diml". Together with (|3.4|) this implies 



To prove part (2) of Theorem ll.il we start by assuming that 
(3.5) {&/,9/}(z )>0, 
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and let G T~Ln be a sequence of vectors such that 

\\(T { f N) - M)(u N )\\ = 0{N~°°) and \\u N \\ = 1. 

If we let u be the distribution on P whose Fourier coefficients are the un, then we 
can rewrite this condition as 

(3.6) (Q-XI)(u)eC°°(P). 

We will now quote Theorem 27.1.11. of (Hormander, Vol IV) asserting that the 
Poisson bracket condition (|3.5|) . translated to the corresponding statement about 
o~q, implies that Q — XI is microlocally subelliptic on the set 

S X0 := {p G Z- 7r(p) = x } C T*P, 

with loss of 1/2 derivatives. Therefore, by (j3.6|) . for all s 6 R, u G ^(°) a ^ ever y 
p G What this means is that for each such p we can write: u = uq + u\ 

where u\ G H^(P) and /? WF(-uo), for all s. Therefore the wave-front set of u 
is disjoint from E XQ . 

3.2. Proof of Theorem ll.2L There are two ingredients in the proof, one is a very 
general localization statement (well-known in the theory of /i-admissible ^DOs), 
and the second Hdrmander's results on microlocal subellipticity that we used in 
the previous section. We begin with the localization result, which is of interest in 
its own right. 

Proposition 3.3. Let be a Berezin-Toeplitz operator with principal sym- 

bol f, and let {un G 7~Cn} be a minimizing sequence of the Rayleigh quotients, 

{~~^ — |ji/<ir~~~~~"~ ' ^ ^ 7~(-n}, where \\un\\ = 1 for all N. Then the distribution 
u := YIn u n G C~ 0O {P) has wave-front set contained in the set of points, p G T*P, 
such that 7t(p) G / -1 (A), and therefore the microsupport of {un} is contained in 

Proof. For simplicity of notation, assume without loss of generality that A = 0. 
Let Q denote the operator on P inducing the and commuting with IT, as 

before. The minimizing sequence is a sequence of eigenstates of the non-negative, 
self-adjoint classical pseudodifferential operator of order zero S = Q*Q on P. We 
will denote by S^ N ' : TCn —> 'Hn the restriction of S to TCn- By assumption un is 
an eigenvector of corresponding to the smallest eigenvalue. Let x G Cq~ 
be a test function, R a zeroth-order ^DO on P, and consider the traces 



Tg = Tr J X (t)Re' itNSN dt. 



If we write the eigenvalues and eigenvectors of SV in the form: 
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where = dim'Hjv an d {ip? }j is an orthonormal basis oITLn, then 

d N 

Just as in the proof of part (0) of Theorem ll.il by taking coherent states as trial 
functions one obtains the estimate: E± = 0(1/N). 

Lemma 3.4. If the microsupport of R is disjoint from the characteristic set ofQ, 
then TjJ2 = 0(N-°°). 

Proof of the Lemma: The operator e~ ltDeS is a Fourier integral operator, and a 
simple wave-front set calculation shows that the wave-front set of the operator 

nos r =no J x (t)e~ itD ° s dt 

is contained in the set 

{(p,p') eZxZ; o-q(p) = = a Q (p')}. 

Therefore, if the microsupport of R is disjoint from the characteristic set of Q, the 
operator R o n o S x is smoothing. Consider next the generating function of the 

T (A0 ; 

T(s) := e lNs Tg2 = tr (u(e is ) oRoUoS x ), 

N 

where U(e ts ) is the operator on P given by composition by the action of e~ ts . 
Another wave- front set calculation shows that T G C co (S 1 ) because Ro II o S x is 
smoothing, and therefore the Fourier coefficients of T are rapidly decreasing. This 
proves the lemma 

Continuing with the proof of the Proposition, choose x so that x > ( we can 
even take it so that x 1S equal to one in a neighborhood of zero), and choose R of 
the form: R = F* o F where F is a zeroth order ^DO on P. Then all terms in the 
sum defining T^l are non-negative, and the previous Lemma implies that, if the 
microsupport of F is disjoint from Char((5), one has: ||-F(uat)|| = 0(N~°°). But 
this implies that, for any such F, 

F(u) = J2HuN)eC°°(P). 

N 

Since we can pick F microlocally elliptic in a neighborhood of any point in the 
complement of Char(Q), the wave-front set of u must be contained in Char(Q). □ 

Turning to the proof of Theorem 11.21 we note that under its assumptions we 
are in a position to apply the subellipticity results of Hormander. Note first that 
any given repeated Poisson bracket of the real and imaginary parts of / evaluated 
at a point x G X is equal to the same repeated Poisson bracket of the real and 
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imaginary parts of the principal symbol of Q, evaluated at any p G Z such that 
7r(p) = x. 

Next, we claim that the hypotheses of Theorem 11.21 imply the hypotheses of 
Theorem 27.1.11. of namely: 

(A) For every x G / _1 (A) there is some j < k and some z G C, 

(3-7) ((H R(z/) )'9(z/)) (*)#() 

where 3^ z ^\ is the Hamilton vector field of $t(zf), considered as a differ- 
ential operator. 

(B) The repeated Poisson bracket above is non-negative (it is positive, actually) 
if j is the smallest integer such that (|3.7|) holds for some z. Moreover such 
j is odd. 

Part (A) follows from hypothesis (2) in Theorem 11.21 bv Corollary 27.2.4 in |12j . 
Lemma 5.1 of (HI shows that hypothesis (1) of Theorem 11.21 implies Hormander's 
Condition (vf), which, as indicated in the remark after Theorem 27.1.11 of |12j . 
implies (B). The proof of this remark is further detailed in the first paragraph in 
the proof of Theorem 4 of [H] . 

Let C Z be the set of points in Z projecting to / _1 (0) (recall that we took 
A = 0). We conclude, by Theorem 27.1.11 of ^2] an d the previous considerations, 
that for all p G O, Q is subelliptic at p with loss of at most 5 = k/ (fc+1) derivatives. 

By Lemma 27.1.5 of |12| . for each p G there exists a zeroth-order \I/DO, A p , 
non-characteristic at p and such that 

Vg G C°°(P) \\A p g\\ { _ 5) < C p (\\Qg\\ L , + ||<7||(-i))- 

An examination of the proof of this lemma shows that we can take the symbol of 
each A p to be non-negative. Since is a cone with compact base, there exists an 
integer K such that the sum of K of the operators A p , call them A\, . . . ,Ak, is 
non-characteristic at each p G 0. We denote such a sum by A = Ylf=i Ar Then 

(3.8) V 5 G C°°(P) \\Ag\\ { _ 5) < C(\\Q 9 \\ L 2 + \\g\\ { _ 1} ) 

for some fixed C > 0. By averaging with respect to the S 1 action on P (which 
preserves Sobolev norms), we can further assume without loss of generality that 
[A,D e ]=0. 

Let B be a microlocal parametrix of A in a neighborhood of such that 
[B, Dq] = 0. Since, by Proposition 13.31 the wave-front set of u is contained in 
0, we have: 

u = BA(u)+g 7 with g G C°°(P). 
By taking Fourier components it follows that 

VN \\u N \\(..g) < C b \\A(un)\\(-S) + \\9n\\(-5), 
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where Cb is the -£^(-<$) norm of B. This together with (|3.8|) implies that for all A~ 

\\Q(un)\\l2 > Ci\\un\\(_ S ) - C2\\gN\\(-6) - ||ujv||(-1). 

Since g £ C°°(P), \\gN\\(-S) = 0(N^°°), while equation (|2"3|) gives us that 
Ci|hiv||(-5) - = ^(N 2 + N + \)~ 6 / 2 -(N 2 + N+ I)- 1 / 2 

Since 5 = k/(l + k) < 1, this proves Theorem 11.21 

4. Examples 

We now look at specific examples of quantized Kahler manifolds, X, and of 
Toeplitz operators. The corresponding Hilbert spaces have canonical bases and 
therefore the Toeplitz operators become sequences of matrices of a specific type 
that we compute. 

4.1. A preliminary example. We begin with a concrete example associated with 
X = C, the plane with its usual complex structure. Although this X is not compact 
(and therefore it does not fit precisely into the general framework of this paper) we 
will "cut it" to the unit disk, both symplectically and quantum-mechanically. This 
leads to a sequence of large matrices to which microlocal methods apply, provided 
one stays away from the boundary of the unit disk. We only consider an explicit 
operator which is the microlocal model of the general case. It will be clear that 
what we do easily generalizes to other operators in this setting. In this section we 
want to be explicit and avoid using the general machinery. 

Recall that the Kahler quantization of the plane gives rise to the Bargmann 
spaces 

Bn = {/:C^C entire ; ||/|| 2 := - f \f(z)\ 2 e~ Nz ~ z dxdy < oo} 

7T Jx 

where z = x + iy and > 0. Elements of arise from the general Kahler 
quantization scheme applied to C. The quantizing line bundle L — > C is holomor- 
phically trivial, and so its sections can be identified with entire functions on C. 
The Hermitian structure on L, however, is not trivial. We introduce the following 
notation for the length function of tp G Bn as a section of L N — > C: 

(4.1) \^{z)\ s :=\^{z)\e- Nz ~ z l 2 . 

Notice that then the norm of ip is the integral of the function \ip\ s with respect to 
the area form. 

A fundamental operator on Bn is the harmonic oscillator (shifted by 1/2, for 
convenience) , 

Op(flH.o.) = N~ x z^- 
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which is a Berezin-Toeplitz operator with symbol H B O (z, z) = zz. The eigenfunc- 
tions and eigenvalues of Op(iJ H .o.) are: 



so that 



N~ 1 z^-z j = N^jzi, and \\z j II = AT- 1 0'+ 1 )/ 2 ^1, 

/v(fc+l)/2 , 

\k) = —j=-z\ fc = 0,l,... 



is an orthonormal basis of £>at. The unit disk is the region of phase space where the 
classical energy z~z is less than one. The analogous object quantum-mechanically is 
the span of the eigenfunctions with eigenvalue less than one, that is the monomials 
z j with j < N. Thus 

TLn = {polynomials in the complex variable z of degree < N } 

where it is now natural to restrict N to be an integer. This setting is close to the 
case of the sphere to be considered in the next section, provided one does not get 
too close to the boundary of the unit disk. The Hilbert space Wat of the sphere can 
also be identified with the space of polynomials in a complex variable of degree at 
most N, although on the sphere \\zj\\ is essentially (C 3 N )~ l l 2 . Symplectically, the 
sphere is the disk with its boundary collapsed to a point. 

The example we will be studying is based in the following observation. Suppose 
one has a Berezin-Toeplitz operator, Q, and a state ip such that = 0- Suppose 
one has a "good" semi-classical cut-off operator, ®n (a projector), which is semi- 
classically the identity in a certain region of the plane (referred to as the allowed 
region), and let P/v = &nQ&n- Then 

= Q N QQ N (il>) + e N Q@^) + @nQ&nW + 0]vQO]v(V) 

where 0^ = I — @n is the complementary projection. The second and third terms 
in this sum will be very small due to the assumed localization properties of @n- It 
follows that 0jv(VO is a good pseudo-mode for P/v (with pseudo-eigenvalue zero) 
if 



IISatVII 

is small, that is, if tp concentrates in the classically allowed region of 0j\r. We will 
take /v to be the orthogonal projection 

0/v : B n — > 7~Cnj 

for which the classically allowed region is the interior of the unit disk. We now 
proceed to make this statement more precise. Our basic tool is the reproducing 
kernel of Bn, in the form of the coherent states: For each w G C, let 

(4.2) <p w (z) := Ne NzW G B N . 
These states have the reproducing property 

(4.3) WeB N , zeC = (il>,<Pz)- 
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Notice that, in particular 

(4.4) \\ip w \\ 2 = (ip w ,<p w ) = Ne N ^ 2 . 

Lemma 4.1. For every 5 > there exists n > such that for all N > n and for 
each w £ C such that \w\ < 1, we have: 

(4-5) ||e^f < i±i||^|| 2 |^| 2 ^e-^-H 2 )V2. 

V Z7T 



Proof. From the Taylor series expression for (p w and the orthogonality relations of 
the monomials z 3 , one can show that 

(Art llft-L II 2 II II 2 (l ^(N + l,N\w\ 2 ) 

(4.6) iiOjvvwii = wvjW y 1 

where F(n, x) = J" °° t n ~~ 1 e~ t dt is the incomplete gamma function. To estimate the 
quantity in parenthesis, notice that 

i-N\w\ 2 

N\ = T(N + 1) = F(N + 1, N\w\ 2 ) + / e"* dt. 

JO 

Therefore, dividing by A"! and making the change of variables: s = t/N, we see 
that the quantity in parenthesis in (j4.6|) equals 

i\tN+1 Aw\ 2 N N + 1 r\ w \ 2 

A 1 -^ S »ds = ^^e- N / e- N ^ds 



N\ Jo Nl 

where f(s) = s — log(s) — 1. This is a decreasing function on s 6 (0,1), and 
therefore the last integrand is maximal at s = \w\ 2 . On the other hand, it is 
elementary to check that 

Ve€(0,l) f(l- e )>e 2 /2, 



and therefore 



A^ 1 f M _V. v , . 



-i 



e 



^ s Jv ds < ^ e -iV . ,2 e -tf(l-M»)V2. 



AT! 7 ~ iV! 

An application of Stirling's formula finishes the proof. □ 

As we now see this Lemma implies the localization properties of Oat: 
Corollary 4.2. For each ip £ B]y, and for each z such that \z\ < 1, 

(4.7) < WW 1*1 N3/ ' e-^H-lW. 
Therefore, for all e > i/iere exist C, a > swc/i i/iai /or all ip £ Bn 

(4.8) //" |e^)(*)l2 dxciy < C HVH 2 iV 3/2 e- aJV . 

.A/|z|<l— e 
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Proof. By the reproducing property (|4.3|) . we have: 

and therefore, by the Cauchy-Swartz inequality and (|4.5[) . for all sufficiently large 
N 

Using (|4.4j) we obtain ()4.7jl . and 1)4. 8 J) follows by integration. □ 

We will need one more general fact about the projection, @n' 
Lemma 4.3. For all z £ C and tp S T~Cn, one has: 

(4.9) |(e^)(z)| a <(iV + l) max |V(e^)| s . 

Proof. This is a consequence of the following formula for Qn- 

(4.10) (e^J (z) = T J2 / e ' lkt ^ tz ) dt > 

together with the fact that the Hermitian weight function e ~ Nzz / 2 [ s S 1 invariant. 

□ 

We now show how the above results can be applied to to a concrete example, 
the sequence P/v = QnQQn where Q is the model operator: 

^ 1 d 

q= nTz + » z 

acting on the space TLn- The symbol of {-P/v} is the function / = [iz + z restricted 
to the unit disk. Notice that {3?/ , 9/} is identically equal to fi 2 — 1, which we 
will assume is negative, i.e. we now take \i 6 (—1, 1)- Notice that for each A £ C 
there exists a unique zq such that A = f(zo). 

Every complex number is an eigenvalue of the operator Q. Specifically, for each 
zq £ C consider the state 

^ zo (z) := ViVV^I 2 / 2 e ™ e -f^-^ 2 . 
One can verify that tp^zo e &n because \y\ < 1. The state (f^ Zo is the (quantum) 

t JV 2 

translate of the basic "squeezed state" at the origin, y Ne~^^ z , to the point zq, 
and its norm is a universal constant (independent of N and of Zq). It is trivial to 
verify that 

Q<p», Z0 (z) = f(z )(fi^,z (z). 
As a section of L® N — > C, the length of (p^ Zo at z is 

(4-11) \^ Z0 (z)\ s = VNe- NF ^ 

where F is the real quadratic function 

F{z) = -&(zzq - -\i{z - z ) 2 ) + z-z/2 + |z | 2 /2- 
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F is non-negative, vanishing exactly at z = zq. The states <f^ Zo are the model of 
the polarized Hermite states associated with a point. Our result is as follows: 

Proposition 4.4. With the previous notation and if \i G (—1,1) and \z$\ < 1, 

|| (P N - \)O N (tp^ Z0 )\\ = Q, N l/4 e -aNj 

||®iv(</V,2o)ll 

for some a > 0. 

Proof. We first will show that 

(4-12) ||©aK^o)I| 2 = ||^J 2 + 0( e -^). 

As noted above, ||</?/z,z || is a constant independent of N. For simplicity, a will 
denote a positive constant that may not be the same at each occurrence. 

Let A denote a disk of radius less than one containing zq in its interior. Then, 
by (|4.9|) and the decay properties of <Pp t z 

(4.13) ||6iv(¥Wo)ll 2 = fj 18^(^)12 dxdy + 0( e - aN ). 

Next, notice that on A, \ip^. Zo (z)\ s is uniformly bounded by a constant times y/~N, 
and therefore (by ()4.9|) again) \QN((p flyZ0 )(z)\ s is bounded by a constant times iV 3 / 2 
there. It follows that 

\®n(^ Z0 )(z)\ 2 s - \^, Z0 (z)\ 2 s = \&h(^, Z0 )(z)\ s ■ 0(N 3 / 2 ), 
with a constant uniformly on z £ A. Therefore, by (|4.8I) . 

@N{tPtM,z )\ 2 sdxdy = // ((^J 2 dxdy + 0(e~ aN ). 
A JJa 

Once again, by the localization properties of l -p lJ , )Zo we have 

2 _ ff I,. |2 > , . , ^t-aN\ 



\Wn,za\\ = JJJ^f,,z \ s dxdy + O{e' 

The last two equations and (|4.13|) imply (|4. 12|) . 

Let us now turn our attention to the vector (P]y — X)Qjy(ip^ tZ0 ). In the standard 
orthonormal basis {|fc}}, the matrix T = (i/ m ) of Pjy is tri-diagonal. Specifically, 
the only non-zero elements of this matrix are: 



(4.14) t k<k+1 = y/k/N and t k+lyk = fiy/k/N . 

Let ^,20 = J2kLo a k\k)- Since ^>^ ZQ is an eigenfunction of Q with eigenvalue A, 

(P N - X)Q N ((p^ ZQ ) = (/iajv-i - Xa N ^j\N) 

We claim that both ajv_i and a at are exponentially small in N. Indeed 



1 



2tt 



a N \N) = — I e- liyt ^ Z0 (e lt z)dt, 
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which implies that for all z S C 



\a N z | < N{N+ i )/2 ™fX-\<Pw>K ze ■ 

where the absolute value is the standard one. Evaluating both sides at z = 1 and 
applying Stirling's formula we obtain 

M < CW-V4 e -iv/2 m a X |^ 2o( ^)|. 

But e _Ar / 2 max i |</^, Zo (e**)| = max 4 |w,«o( ei *)ls> an d therefore 
e-^maxl^Ce*)! = 0{VNe~ aN ) 

since |zo| < 1 fsee l4.1ll and the remarks following it). Therefore = 0(N 1 l i e~ aN ), 
and similarly for ojv— i- D 

Much more generally we can start with a pseudo-mode of a general Berezin- 
Toeplitz operator, Q, on Bargmann space. By the localization properties of the 
projector Oat, the projection by ®n of the pseudomode will be a pseudomode of 

e N Qe N . 

4.2. Quantization of the torus. We consider the standard torus, X = C/A, 
A = Z 2 , with the complex structure arising from that of C. A quantizing line 
bundle on X is holomorphically trivial when pulled-back to C, therefore its sections 
can be identified with entire functions on C satisfying a transformation law with 
respect to translations by elements of A. It is well-known that the functions that 
arise in this manner are precisely theta functions. The quantizing line bundle is 
not unique, since one can always tensor a given one with the flat line bundles over 
X. This gives rise to theta functions with characteristics. 

4.2.1. The Hilbert spaces. A quantizing line bundle over X can be constructed 
from a cocycle x '■ C x A — > C \ {0} given by: 

(4.15) X (z,m + in) = (-\) mn e 'r[«(m-m)+i(m 2 +n 2 )] e -27ri[m M + W ] ) 

where [i and v are fixed real numbers (the so-called characteristics of the bundle). 
X is called a cocycle because it satisfies the condition 

(4.16) x(*. A) x (z + A, ft) = X (z, A + /x)- 

The quantizing line bundle is the quotient of C x C by the equivalence relation: 

(z, a) ~ (w, b) 44> 3A G A such that (w, b) = [z + A , x(z, A)a). 

For simplicity we will only consider here theta functions with characteristics (fj,, v) = 
(0,0). _ 

We observe the following features of this construction: 
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(1) The sections of this line bundle, L, are naturally identified with the func- 
tions / : C —* C such that 

(4.17) V(z,A)£CxA f(z + X) = X (z,X)f(z). 
(Indeed the section associated to one such / is defined by: 

•/(M) = [(*,/(*))] 

where the square brackets denote equivalence classes.) 

(2) For any integer N the iV-th power of x, \ N ls again a cocycle. The line 
bundle it defines is the iV-th tensor power of L, L® N . 

(3) A Hermitian structure on L is defined by a function h : C — > M + satisfying: 

(4.18) h(z) = \ x (z,\)\h(z + X). 

The Hermitian metric we will consider is: |[(z,a)]| = |a|e _7r l 2 l 2 / 2 . 

Definition 4.5. The space Wat of holomorphic sections of the line bundle L® N is 
the space of entire functions / : C — > C satisfying: Vz G C, m + in G A 

/(z + m + in) = (_i)^" e ^W™-»)+i(" l2 +« 3 )]/( z ), 

Its Hilbert space structure is given by the inner product 

(4.19) <f,9>= I f(z)^)e- N ^ 2 dxdy, 



IT 

where T is a fundamental domain for A. 

The transformation law (|4.17jl is not the standard one for theta functions (see 
0, d). However, if / € H N , then 

f x (z) := e~ N ™ 2 / 2 f(z) 
satisfies the classical transformation law 

f x (z + m + in) = e NAn 2 -2inz) f x( z y 

For future reference we also associate to / G Ti^ the function 

f y (z) :=e N ™ 2 / 2 f(z) 

which satisfies the transformation law 

f y (z + m + in) = e N *( m ' 2 + 2mz ) fy( z ). 

Note that, in particular 

r(z + n)=f x (z), and f{z + im) = f y (z). 

Exploiting these periodicity conditions we now exhibit two (dual) basis of Wat. 
We begin with the functions f x . They can be expanded in Fourier series, 

f X ( z )= ^2 a me(mz) 
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where we let e(z) := e 2mz . The transformation law for f x becomes a relation 
among the Fourier coefficients, namely a m +Ar n = e _7r ( Arn + 2mn ) a m . This shows 
that the dimension of the space of theta functions of order N is N, as the values 
of ao, • • • o/v-l determine the Fourier series. This leads to considering some special 
theta functions obtained by letting exactly one of the coefficients ao, . . . ajv-i be 
non-zero. These functions give rise to a basis of Wat. More precisely: 

Lemma 4.6. For j = 0, . . . , N — 1, let m N \z) be defined by the Fourier series 



j 

oo 



$f\z) = (2iV) 1 /4 e A f - 2 /2 £ e-^ n+ i^ 2 e(z(j + Nn)). 

n=— oo 

Then £ Wat, and f/ie se£ , i = . . . , N — 1} is an orthonormal basis of 
Hn- 

We can carry out a similar construction by considering the Fourier series of the 
functions f y where / £ Hn- This results in a different basis of Hn- 

Lemma 4.7. For k = 0, . . . , N — 1, let p[ N \z) be defined by the Fourier series 



ff\z) = (2iV) 1 / 4 e-^ 2 / 2 e- nN{m+k/N)2 e(-iz(k + Nm)). 



'k 

oo 



m=— oo 



Then /^i £ Wat, and i/ie sei , = . . . , N — 1} is an orthonormal basis of 
Hn- 

It's a beautiful fact that the matrix relating these two bases is the discrete 
Fourier transform. We claim that 

JV-l 

(4.20) { k N \z) = Y, e~ 2 ^/ N $f\z). 

3=0 

We refer to and especially |15j . Proposition 3.17 for the reason for this. 

4.2.2. Matrix coefficients. We now turn to a calculation of matrix coefficients of 
Toeplitz operators on the torus. We recall the following result of (Corollary 
4.8): 

Lemma 4.8. Let (x, y) be standard coordinates on the torus, so that z = x + iy 
in the previous formulas. Let f(x,y) be a symbol that is in fact a smooth periodic 
function of y alone. Then the #j are eigenvectors of the B-T operator T^ , 

oo 

T?$f) =X fhf\ where xf > = £ a n e~™ 2 1™ e~^/ N ~ f\ y= _ j/N 

n=— oo 

and where the a n are the Fourier coefficients of f (with respect to y). 
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Notice that in particular all such operators are normal. Similarly, the basis 
consists of eigenvectors of any Toeplitz operator with total symbol a function of x 
alone. 

Continuing our calculations of matrix coefficients, let us now take a symbol / 
of the form: 

(4.21) f{x,y) = h(y)e 27rilx , where I G Z. 

Lemma 4.9. Iff is of the form \4-21\) , then (f'&j , ) is zero unless k = (j + l) 
mod N, in which case 

More precisely, the matrix coefficient above is equal to 

g — 7rJ 2 /27V {e -*Ay/2N h){ _2j+l : 



where A y = ^j. 

Proof. It is easy to verify the first statement, and, in case k = j + I mod N, one 
computes that 

where 

oo 

tfj-f(y) = e _7r ' 2 / 2Ar ^ e -27r7V[y+n+(2j+0/2iV] 2 _ 
n=— oo 

If /i = 1, then we have: 

( e 2^^W^W )roodiv)= ^]V e -^/2iV £ / e -2^[, + n + (2 J+ 0/2^V] 2 dy 

,V1 J 



/OO 
e -2-^ 2 ds = e~^/ 2N . 
-oo 

This is 1 + 0(1/N) and therefore satisfies the desired estimate. To proceed in 
case h is not constant, notice that by the Poisson summation formula Vl/j i can be 
written as 

c -irl 2 /2N 00 



?2N , 

k=—oo 



Therefore 



k=—oo 

□ 
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Corollary 4.10. Let k\ < &2 be two integers, and let hi(x) be smooth 1-periodic 
functions, ki < / < k-i- Then there exists a B-T operator on the torus, {T^}, 
with principal symbol 

l=k 2 

f(x,y) = ^hi(y)e 2mlx 
l=ki 

and such that for each N the matrix entries ofT^ in the basis {$^} are given 
by the formula 

(4.23) (TW$f\# ij+l)modN ) = M"^) + 0(N~™), 

where the estimate is uniform (in j, I). 

Proof. For each I there exists a sequence of periodic functions fi >m , m = 0, 1, . . . 
such that 

oo 
m=0 

where the left-hand side is considered a formal power series of 1/N (clearly fio(y) = 
hi(y)). By the Borel summation method, there exists an iV-dependent function 
fi(y,N) periodic and smooth in y such that fi(y,N) ~ Ylm=o N ~ m fl,m{y), esti- 
mates the C°° topology. By the previous Lemma, the matrix coefficients of the 
Toeplitz operator with iV-dependent multiplier f(x,y;N) = Y^,i=k fl G/> e 2 ™ 1 ® 
satisfy (|4.23j) . Since the sum over I is finite the estimates are uniform. □ 

4.2.3. Localization of pseudomodes. We now proceed to describe in concrete terms 
the phase space localization of the pseudomodes constructed in §2. We begin 
by pointing out that the concentrate, as N — ► oo, along the circles on the 
torus defined by y =constant. This is because the flj are eigenvectors of Toeplitz 
operators with symbols / = f(y). A precise statement is: 

Lemma 4.11. Fix yo = j^- a rational number modulo 1. Then the microsupport 
of the sequence {fi^ ^ ', fe = 1, 2, ...} equals the circle y = yo on the torus. The 
j3j accumulate along circles x = xq in an analogous fashion. 

An even more precise statement is that the sequences of this lemma are Legen- 
drian states associated to the corresponding circles, in the sense of jB]. 

Consider now a pseudomode, f/w, associated with a Toeplitz operator Ti N ' on 
the torus, with microsupport (xQ,yo) £ X. It follows from the previous lemma 
that the components of i^n 

in the basis of the concentrate around the values j where j/N ~ yo. Similarly, 
the coefficients 
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concentrate around the values j where j/N w x . More precisely, the concentra- 
tion occurs in a neighborhood of size 0(s/~N) of these values. As mentioned, the 
sequence is the finite Fourier transform of Oj N ^ . Therefore, in these coordinates 
the pseudomode is localized on both sides of the Fourier transform. 

4.3. The complex projective line. We begin by reviewing the quantization 
of P 1 , for completeness. The quantization of the complex projective line arises 
in connection with the irreducible representations of SU(2). Recall that up to 
isomorphism such irreducible representations are those realized in the spaces 

Hn := { f(wi,W2) ; / a homogeneous polynomial of degree N }. 

Specifically, if / G Hn and g G SU(2), then 

(g ■ f)(w 1 ,w 2 ) = fig' 1 ■ (w 1 ,w 2 )) 

where the action on the right-hand side is the natural action of SU(2) on C 2 . Here 
N = 0, 1, 2, ■ ■ ■ is a non-negative integer. These representations are unitary if we 
put on Hn the Hermitian inner product 

</l,/2>= / hT2dV s3 , 

where S 3 C C 2 is the unit sphere and dV S i is its volume form. We note without 
proof that the vectors 

(4.24) | j, N) = ^^±±^clw{ w *-i, 0<j<N 

form an orthonormal basis of Hn (consisting of eigenvectors for the operator in- 
duced by o"3 G su(2), see below). 

The circle group, S 1 C C, acts freely on S 3 , by complex multiplication. There- 
fore S 3 is a circle bundle over the abstract quotient, X := S 3 /S 1 , which can be 
identified with the space P 1 of all complex lines through the origin in C 2 . There is 
a natural Hermitian line bundle, L*, over P 1 (the one whose fiber at £ G P 1 is £ it- 
self), and in fact P = S 3 C L* is the unit circle bundle. The dual bundle, L —> P 1 , 
is the so-called hyperplane bundle and, as it turns out, quantizes P 1 (meaning that 
the natural connection on it has curvature the Fubini-Study symplectic form of 
P 1 ). By homogeneity, the functions on Hn transform very simply along the orbits 
of S 1 , and the inner product above is invariant under S 1 . In fact we can regard the 
elements of Hn precisely as the holomorphic sections of L® N . We are therefore 
exactly in the setting of the previous sections. 

The space P 1 is isomorphic to a two-dimensional sphere, as follows. Define a 
map: 

(4.25) $ : P 1 su(2) 

by the following rule: For each f £ P 1 , is the matrix having £ as an eigenspace, 
with associated eigenvalue i/2, and £ L as another eigenspace, with associated 
eigenvalue —i/2. (The choice of spectrum is dictated by the normalization that 
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the area of S 3 /S 2 agrees with the one induced by the Killing form.) One can show 
that, if we write £ = [wi,W2] € P 1 , where (w\,W2) G S 3 is a representative, then 
the previous map is: 

It is a fact that $ is a moment map for the natural SU(2) action on P . Recall 
also that the matrices 

1/0 A 1 fO -l\ 1 (i \ 

form a standard orthogonal basis of su(2) (such that [a±, c 2] = 03, etc., and ||<jj|| = 
1 /2) if we give to su(2) the SU-invariant inner product 

(A,B) = -±TrAB. 

Clearly <& is an equivariant diffeomorphism onto its image, which, geometrically, is 
the sphere of radius 1/2 and, algebraically, a (co)adjoint orbit. (It turns out that 
the symplectic form on P 1 is twice the area form. In general the symplectic form 
on an orbit of radius s is the area form divided by s, see ^S] Pg- 54.) We will 
henceforth identify P 1 with this sphere/co-adjoint orbit. 

Let Xj : su(2) — > R be the j-th coordinate function, Xj(A) = (A,2aj), so that 
the image of $ is the sphere X^i=i x | = \- ^ rie description of P 1 as a sphere means 
that we can speak of restrictions of linear functions from su(2) to P 1 . In particular, 
a crucial role in what follows will be played by the function I : P 1 — > [0, 1] given 
by 

(4.28) i(t) = (m, ^) + \ 

The Hamilton flow of this function with respect to the natural symplectic structure 
on P 1 (as a co-adjoint orbit) is given by the action of the one-parameter subgroup 
exp(2to"3), and geometrically is rotation around the 03 axis. Accordingly, we in- 
troduce the polar angle, 8, regarded as a multivalued function on P 1 (undefined 
at the poles). (1,0) are action-angle coordinates on P , and, in particular, the 
symplectic form on P 1 is 

lo = dl A d6. 

Any Berezin-Toeplitz operator gives rise to a sequence of matrices, namely, the 
matrices representing the operator in the canonical basis, {\j, N)}, of TLn- We will 
compute below (approximately) the matrix of a Toeplitz operator with principal 
symbol a given function / : P 1 — > C. We begin with some remarks about such 
functions. Let 

00 

(4.29) f(I,6)= ^flW 

i=— 00 
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be the Fourier series expansion of / with respect to the action of S 1 by rotations 
around the 03 axis. Although we are writing this expansion in in action-angle vari- 
ables, each summand is a smooth function on the sphere, which imposes boundary 
conditions on the f\. Specifically, one can prove that the functions fi(I) such that 
fi(I)e lW is smooth as a function on P 1 are of the form 

(4.30) fi{x) = (x(l - x)) 1 / 2 gi(x) 

where gi has a smooth extension to a neighborhood of [0, 1]. In fact, x\ + 1x2 = 

(1(1 -I)) 1 ' 2 e iW = ( Xl + ix 2 ) 1 . 

Therefore, if dEP) holds, 

h(iy w =g l (I)(x 1 +ix 2 ) 1 

which is clearly smooth on the sphere if gi £ C°°[0, 1]. 

Let us now turn to the computation of matrix elements of B-T operators on 
P . By linearity of Toeplitz quantization, it will suffice to compute the matrix of 
a Toeplitz operator with symbol e lW fi(I) for a given integer I. We begin with the 
case / = 0. 

Lemma 4.12. Let a be a smooth function on [0, 1]. Then there is a B-T operator 
on the sphere, {A^}, with principal symbol aol , which is diagonal in the standard 
basis of TCn and whose j-th diagonal entry is a(j/(l + N)), < j < N. 

Proof. The sequence of operators, Z = {Z^ : TCn —> TCn} such that 

zW\j,N) = j^\j,N), 0<j<N 

(where \ j, N) is defined in (|4.24l) is a Berezin- Toeplitz operator with symbol / (see 
Lemma 3.4 and the ensuing discussion in .5!). Let 

2tt J 

where a is the Fourier transform of a compactly-supported smooth extension of 
a. Since Z is a self-adjoint Toeplitz operator of order zero, is a unitary 

Toeplitz operator of order zero (see for example Proposition 12 of [7]) and symbol 
e ltxs . It is easy to check that {A( N '} has the desired properties. □ 

To compute the matrix elements of operators with a smooth symbol e tW fi(I) 
where I 7^ 0, we introduce the raising and lowering operators. Let 

:= % x the operator induced by on TCn, k = 1,2, 3. 

Then: 

- 1 / d d \ N if d d 



dw\ dw2 ' 2 V dw\ dw 
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and 



N l( d d 

J 3 = -\Wi- W2- 



2 V dw\ dw2 ■ 
In particular, the lowering and raising operators, J± = ± iJ%, are: 

J+ =Wl- , J_ =w 2 - 



dw2 ' dw\ 
(and the vectors \j), j = 0, . . . ,N are eigenvectors of , with eigenvalue j — ^-^-.) 

Lemma 4.13. The matrix of , resp. , in the standard basis has zero entries 
except along the supra- diagonal, resp. infra- diagonal, along which the entries are 
equal to 



^j(N-j + l), j = l,...,N. 



Moreover, the sequence {jfJ±} is a B-T operator with symbol x\ ± ix 2 - 

The first two statements follow a simple calculation; for the last statement we 
refer to 

We are now in a position to describe the matrices of the B-T operators on the 
sphere: 

Proposition 4.14. Let f : P 1 — > C be a smooth function with a finite Fourier 
series, {^.29j) , in action-angle variables, and let fi{x) = (x(l — x)) 1 ' 2 gi{x) with 
gi £ C°°[0, 1]. Let be the matrix of J± , described in the previous Lemma. 

Then there exists a B-T operator, {T™}, with symbol f and such that the sequence 
of matrices {T^} of {T^} in the standard basis ofTL^ satisfies: 

(4.31) TW = Y.{ M % ] n^ A{N) ^)+0(N-n 

l 

where A^ N \gi) is the (N + 1) x (JV + 1) diagonal matrix with diagonal entries 
gi(j/(N + 1)) and the estimate is in any matrix norm. 

Proof. By linearity and the assumption that the sum (|4,29|) is finite it suffices to 
prove the proposition for / of the form /(J, 0) = e lW fi(L). The case I = is 
covered by Lemma 14.121 It suffices to consider the case / > 0. Applying Lemma 
14.121 to gi we obtain a diagonal B-T operator, {^4^ (<#)}, with diagonal entries 
gi{(j - 1)/N), < j < N. It is clear that the B-T operator 

1 JN\\ l 



T (N) := (±jwy oA (N) 



has the desired properties. Notice that one can choose any order of the products 
appearing in 1)4.31(1 and still find a B-T operator with the desired properties. □ 
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V = exp(ito"2) 



4.3.1. Linear Hamiltonians. Recall the moment map, <E> : P 1 — > su(2) given by 
(|4.26|) . Given M E sl(2,C), we can pull-back by the complex- valued linear 
function on su(2), 

su(2) B -~Tr(AM). 

Let us denote the pull-back by Fm '■ P 1 —>■ C; specifically 

W G P 1 F M (£) = -- Tr($(£)M). 

We will call functions such as Fm linear Hamiltonians. Since $ is a moment map, 
the assignment M — > Fm is a Lie algebra morphism: 

VMi, M 2 E sl(2, C) F [Mu m 2 ] = {F Ml , F M2 }. 

In particular, if we continue to denote by Xj : P 1 — > M the restriction to P 1 of the 
coordinate functions, then we have the identity: {x\ , 22} = %3, and also its cyclic 
permutations. 

If M is semi-simple, there exists g E SU(2) such that gMg~ 1 is diagonal, i.e. 
3g G SU(2), M G C \ gMg^ 1 = n a 3 . 
Here's a very concrete example. Take 

cosh(t/2) -isinh(i/2)\ 
v isinh(t/2) cosh(i/2) J 

and consider A = A{t) G sl(2, C) equal to 

A = VaiV^ 1 = i sinh(t) a\ + cosh(t) 03. 

The classical Hamiltonian, Fa '■ P 1 —>■ C is i*A = isinh(t)xi + cosh(i) x%, and its 
image is the interior of an ellipse, 

(a 321 — I y~ < - 

{ ' cosh(t) 2 sinh(t) 2 - 4" 

If we let T^ N ' : Ti^ — > H-n be 1/N times the operator image of jA by the represen- 
tation pn, then, for all t, T^ N > is diagonalizable with real spectrum {jj- — ; j = 
0, . . . , N — 1 }, which is contained in the major axis of the image of Fa- 

Notice that {^R-Fa , 9 -Fa} = cosh(t) sinh(t) X2, and therefore for every point 
in the image of Fa there is exactly one x where this Poisson bracket is negative. 
Therefore, for each A in the interior of the image the norm of the resolvent \\(TW- 
A/)- 1 !! is OiN 00 ). 

The level sets of the norm of the resolvent for this example resemble the equipo- 
tential curves of a uniform electric charge distribution on the line segment [—5, 5]. 
Notice, incidentally, that the image under Fa of the level curves of the Poisson 
bracket {?R.Fa , are ellipses crossing this line segment. This example there- 

fore shows that there is not a direct relationship between the norm of the resolvent 
and the size of the Poisson bracket of the real and imaginary parts of the symbol. 
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Finally, notice that every A on the boundary of the elliptical region (J4.32|) satis- 
fies the hypotheses of Theorem 11.21 The set / _1 (A) consists of exactly one point, 
where {MFa , ^Fg} = cosh(t) sinh(i)x2 = 0. However, one of the double brackets 
involving real and imaginary parts of Fa is non-zero at this point. Thus Theorem 
11.21 applies, with k = 2. Notice that, in spite of the estimates ()1.9|) the distance 
from A to spectrum is 0(1). This is in agreement with Theorem 3 of 0: Thinking 
of P 1 as a real manifold O, the symbol / has an obvious holomorphic extension to 
a complexification of O, namely, the (co)adjoint orbit of SL(2, C) through 03. 

5. Final Remarks 

5.1. On the numerical range. Let Tf = {T^} be a Berezin-Toeplitz operator 
with symbol / : X — ► C For each N, the numerical range of is the set 

(5.1) W N := { (T^ifi, if,) ; V E H N , ||V|| = 1 }• 

We define Woo as the limit of the ranges W/v as N — » 00: 

Definition 5.1. A complex number A is in Woo iff for all e > there exists K > 
such that for all N > K 

A e (A) n W/v 7^ 0, 
where A e (A) is the disc of radius e centered at A. 

Proposition 5.2. Woo is the convex hull of the image of the classical symbol, f . 

Proof. It is well-known that, for each N, W/v is convex (see jllj ). It follows easily 
that Woo is convex as well. Moreover, if x G X and is a coherent state at x, 
then 

This shows that Woo contains the image of /. 

To show that Woo is actually the convex hull, consider a line of equation ax+by = 
c, and let A S W^. By definition, there exists a sequence A at S Wjy converging to 
A, and therefore there exists a sequence of unit vectors {V>tv G T~(-n} such that 

X N = {T^ N ,ip N ) - A. 

Let us write A at = xjy+iyN for the real and imaginary parts of Atv, and / = p±+ip2 
for the real and imaginary parts of the symbol /. Then 

ax N + by N = \ {api + bp 2 ) \iPn\1 dm + 0(1/N). 
Jx 

Assume that the region: ax + by > c does not intersect the image of /. Then for 
all E X a(j)(m) + bg(m) < c, and therefore 

ax N + by N < c / |^at| 2 dm N + 0(1/N) = c + 0(1/N). 
Jx 
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Letting N — > oo, we obtain that ax,^ + by^ < c, where A = + iy oo- Thus A, 
and therefore all of , is on the same side of the line ax + by = c as the image 
of/. □ 

5.2. The weak Szego limit theorem. In the non-selfadjoint case, one has the 
following version of the Szego limit theorem for B-T operators: 

Proposition 5.3. Let Tf = |T^} be a B-T operator with principal symbol f : 
X — > C, and let F(z) be a function of a complex variable analytic on a simply- 
connected region containing the image of f. Then 

where dm is the Liouville measure of X. 

Proof. The idea of the proof is standard; we include some details for completeness 
and to verify that the usual proof is valid in the current setting. As in the theory 
of pseudodifferential operators, for each A not in the image of / one can construct 
a B-T operator, B\, such that the Schwartz kernel of 

:= flf) o (TW - A) - I Hn 

is a smooth section of the bundle Hom(L Ar , L N ) — > X which is rapidly decreasing in 
N (together with all its derivatives). The principal symbol of {B^} is (/ — A) -1 . 
Multiplying on the right by (T^ — A) -1 we obtain 

(TW-xr^Bx + S™ 

where {S 1 ^} has the same properties as ■ Let T be a simple closed curve, 
positively oriented, contained in a region where F(z) is analytic and containing 
the image of /. Then 

The trace of the second term on the right-hand side is 0(N~°°), while 



(N) 



d\. 



dimHN It dim7i/v Jv 

But it is known that Tr b[ N) = ^ j x (f - A)" 1 dm + 0(l/N), where 

the estimate is uniform for A on compact sets away from the image of /. □ 

In particular, if A is a complex number away from the image of /, one has: 

< 5 ' 3 > ^^"-^-WX J x y- x im + W N). 

Clearly the left-hand side of this equation is a sequence of analytic functions in A 
defined away from the union of the spectra of the T^ . On the other hand, the 
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integral J x dm is analytic away from the image of /. Examples show that the 

spectral radius of has a limit, R, such that the image of / is not contained in 
the circle of radius R. It is not immediate to extend Q5.3I) to A with |A| > R but 
inside the image of /. 

Appendix A. Hermite distributions and symbol calculus 

Oscillatory integrals and symbols 

We place ourselves in the setting of §2.2: Let 1Z C Z be a closed conic isotropic 
submanifold, and Y X the reduced isotropic submanifold of X. Hermite dis- 
tributions in I m (P, 1Z) are defined locally as oscillatory integrals. To write down 
an explicit form for these integrals, we'll choose Darboux coordinates (q,p) £ M 2n 
for X with q = (q', q") El'x R n ~ l such that Y = {q" = p = 0}. Here I = dimY. 
For P we then have coordinates z = (q,p,9), and we can always find a function 
h(p, q) such that a\ p= o = d9 — dh\ p= o. The lift of Y to the isotropic 1Z is given by 
specifying that 9 = h(q', 0). 

Now we'll introduce phase coordinates (r, 771,772) £ K+ x M n_ ' x M ra , and a phase 
function 

<p(z, t, 77) = t(9 - h(q,p)) + 771 • q" + r? 2 • p, 
which parametrizes 1Z. A distribution u G I m (P,7Z) can be written locally as 

where the amplitude a(z, r, u) is rapidly decreasing in u and has a expansion in r 
of the form 

00 

a(z,T, V )^T m -y 2 J2 T ' j/2a ^- 
3=0 

The symplectic spinor symbol o~(u) should be thought of as ao written in a suit- 
ably invariant way. The choice of <f> defines for each p E 1Z a canonical isomorphism 
between T*(M. n ~ l x W n ) and the symplectic normal N p , by which the H 00 (N p ) por- 
tion of the symbol may be pulled back to a rapidly decreasing function of m This 
gives the ry dependence of ao, while the half- form portion of the symbol encodes 
the dependence on z € P in a coordinate-independent way. 

The construction of u in fci.ll gives u such that a u = v u ® k u (g> e. Here e is a 
Gaussian in iJ OQ (^°) which can be explicitly computed in terms of the metric. The 
component k u £ H 00 (E p ) satisfies C{n u ) = according to the construction. This 
does not fix k u , but we are free to assume that k u is also Gaussian. Here C is the 
operator on H^Ep) given by the Heisenberg representation of the Hamiltonian 
vector field £ of oq. But according to Theorem 11.4 of 0], £z is just the lift of 
the Hamiltonian vector field Hj on X up to Z. So C could be written explicitly in 
terms of /. 
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By this construction of u, the leading amplitude ao(z,r]) is a Gaussian function 
of rj: 

a (z, V )=c(z)e-^ A ^ 2 

where A is a symmetric matrix with positive definite real part. Arguing as in the 
proof of Theorem 3.12 in Hj. we can write um = Hn(u) as 

u N ((q,p,h(q,p)) =e iNh{q ' p) j e~ im e i<f> a{z , r, rj/y/T) dO dr dr,, 

(near q" = p = 0) and apply stationary phase to the N — * oo limit. The asymptotic 
result is that 

U N (q, P ,h(q,p)) ~ c ( q >) N n-l/2-l/2 e -N( q »,p)A( g »,py/2^ 

where c(q') is independent of N. In particular, at a point z$ G P lying above F, 
we have 

(A.l) u^o) ~ c(z )iV n - i/2 - 1/2 . 

A review of the Hermite calculus 

The composition of Hermite distributions is described in Theorem 9.4 of (I], 
which implies that 

II: I m (R,K) — I m (R,K). 

The symbol calculus corresponding to this composition is based on symplectic 
linear algebra found in §6 of jlj, where the reader can find full details. Since the 
symbol calculus is somewhat involved, we begin with a review of the general details 
of the symbol map before applying it to our case. 

Let V and W be symplectic vector spaces, T C V x W~ a Lagrangian subspace 
and S C W an isotropic subspace. We think of T as a canonical relation from W 
to V; To £ is an isotropic subspace of V. We will make the simplifying assumption 
that 

(A.2) f/ ={w£E;(0,w)er} = 0, 

valid in the applications of the calculus to this paper. 

We assume given a symplectic spinor on S and a half- form on T. Recall that 
if Hoeiy) denotes the space of C°° vectors in the metaplectic representation of 
the metaplectic group of the symplectic vector space V, the space of symplectic 
spinors on S is i^E /!]) A V2 ( S )- 

Under the assumption QA.2JI . the (linear) symbol map of the Hermite calculus 
is a linear map 

(A.3) floo(E7£) A 1/2 (s) » A 1/2 (r) -» #oc((r o s)°/r o s) ® A 1/2 (r ° s). 

Our first goal here is to describe the map (lA~3l) . There are two ingredients in its 
construction, which will be examined separately. First however we must introduce 
the following vector spaces: 
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Definition A.l. U x := { w G E° ; (0, w) G T } C^.and 

U := image of U x in E°/E XJ X . 
These spaces enter the calculus in the following way: 
Lemma A. 2. The subspace 

U C E°/E 

is isotropic, and there is a natural identification 

(A.4) u°/u = (r o E)°/r o s . 

The first ingredient in the symbol map is a canonical isomorphism: 
Lemma A. 3. Under lA.Sfy . there exists a canonical isomorphism 

A 1/2 (E)®A 1/2 (r) = A 1/2 (^n® A 1/2 (r°£). 

Proof. Let 

(A.5) p : r e S -» i/f 

be the map = w; — wi. One can show that the image of this map 

is exactly U°. Moreover, because of (|A.2j) . the projection ((v,w),wi) i— > v is an 
isomorphism 

kex(p) = T o S. 
This is the non-trivial vertical arrow in the diagram: 



I 

-> ker(/)) -> T E ?7f ^ 
(A.6) I 

roE 

I 
o 

The horizontal sequence is just the natural short exact sequence associated to the 
surjection p. 

Having established the existence of these exact sequences the desired isomor- 

1 /2 

phism follows from the behavior of the functor f\ 1 when applied to short exact 
sequences and to direct sums. □ 

The second ingredient in the Hermite calculus is the following: 
Lemma A.4. Under the assumption HA.fy) . there exists a canonical map 

S(E°/£) - A _1/2 (^i) ® S((T o E)°/r o E) . 
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Proof. This is based on Lemma (|A.2[) and the following generalization of a map 
defined by Kostant: 

Claim: Let A (in our case we will take A = be a symplectic vector space, 

and U C A an isotropic subspace. Then there is a natural map 

(A.7) H OQ (A) -> 7\~ 1/2 (U) ® H^/U) . 

The desired map follows from these two claims, if we recall that U = U\ (because 
of (IA~2|) . □ 

To obtain the symbol map (|A..3|) . tensor the maps from the lemmas and use the 
fact that the symplectic form on W defines a natural identification 

A~ V Vi)®A 1/a (Oi) = C- 

Proof of Proposition \2. 51 

Given v G I m (P,7Z), we want to calculate the symbol of u = U(v) G I^(R,TZ). 
In order to apply the symbol calculus reviewed above, we need to rewrite H(v) as 
the composition of a Lagrangian distribution with a Hermite. Thus we introduce 
7r : P x P — > P, the projection through the left factor, and F: PxP^PxPxP 
the map P(f>i,P2) = (pi,P2,P2)- We can then write 

n( w ) = ^*(nsn), 

where (abusing notation slightly) LT here denotes the integral kernel. 

For the symbol calculation it suffices to localize to p G 1Z C T*P. To simplify 
notation we will introduce vector spaces 

V = T p (T*P), W = VxVxV. 

In W two copies of V carry the opposite symplectic form, but to simplify notation 
we just denote the vector space. We define the vector spaces R p , Z p , E p , etc. as 
in H3.ll Note that Z p is a symplectic subspace of V, and R p is isotropic in Z p . 
Locally, the symbol of v G I(P, R) can be written 

<r(v)\ p = v (8) k (8) A 

G A V2 (^p) ® ^oo(^p) ® Hoo(^). 

As an integral kernel, II G /(P x P,Zp^) with symbol 

<7(II) = V^®e®ee A V2 (^p) 8> H^Z;) ® iJoo(^°) ; 

where is the canonical volume form given by the symplectic form on Z p . 
The operator ir*F* is a Lagrangian FIO with canonical relation 

T = {(v; v, w, w); v, w G V} C V x W. 

(To see that it's Lagrangian one needs to keep track of the signs of the symplectic 
forms). The combination UMu is a Hermite distribution associated to the isotropic 

£ = {(z,z,y); z G Z p ,y G P p }. 
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Note that r o £ = R p . 

First we'll describe the spaces that play a role in Lemmas Al-4. To begin, note 
that R° was calculated in Lemma 2.3 to be Z° ffi E p , where E p is the symplectic 
normal of R p as a subspace of Z p . Then we see that 

S°/S = Z p ° xZ° p x (Z° p ®E p ). 

Then 

U! = {(0,v,v); veZ° p }GW, 

and U is the same set as a subspace of S°/E, so that L7 /?7 = Z° ffi £ p = i?°/-R p . 
Introducing the map p as above, we see that 

Image(p) = U{ = XJ X ffi (V x Z p x Z p ) C W, 

and ker p = R p . 

With these identifications, we can decompose the pieces in the symbol map 
(A. 2). First of all, in Lemma A. 3 the map reduces to the identity map on 

A 1/2 (z p )©A 1/2 W©A 1/2 (n 

combined with the obvious decomposition 

A 1/2 (^) = A 1/2 (^)®A 1/2 (^°). 

Also note that since = Z p x Z p x (Z p © E p ) we have 

ffoo(S7S) = #oc(^°) © #oc(^°) © #oc(^°) © ^oo(^p). 

This means we can break the symbol map (A. 2) into three pieces. The first is 
the identity map 

#oo(£ P ) © f\ 1/2 R p - H^Ep) © f\ 1/2 R p . 

The second is a canonical pairing (A. 6) applied to the symplectic space Z° x Z° 
with the diagonal as Lagrangian subspace: 

Hoo(Z°) © Hoo{Z° p ) A _1/2 (^p)- 
Finally, the third is the natural identification 

A~ 1/2 (^°) © A 1/2 (^°) ® A 1/2 (^) ^ A 1/2 (^) = c, 



defined by the symplectic forms. 

Applying this to the symbol of H(u), the third piece shows some natural half- 
forms canceling, the second gives contributes a factor (e, A) which can be absorbed 
by changing the half-form component v ^ v 1 . And the first map then gives the 
stated conclusion: 

a(U(v)) = v © k© e. 
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